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SYMBOLS 


Cf local skin-friction coefficient 

e turbulent mixing energy (two-equation model) 

f ratio of rates of strain in distorted homogeneous flows (eq. (53)) 
h mass-averaged static enthalpy 
h" mass-averaged fluctuating static enthalpy 

K distortion parameter (eq. (63)) 

, ku T 

k+ surface roughness parameter, —j- 

£ length scale (eq. (15)) 

p mean pressure 

p’ fluctuating pressure 

? e production term in turbulence kinetic energy equation (eq. (1)) 

qj mass-averaged turbulent heat-flux vector (eqs . (9)) 

q.T mass-averaged total heat flux vector, molecular plus turbulent 
J mechanisms (eq. (7)) 

R e scaling Reynolds number for A* (eqs. (16) or (39)) 

R^ scaling Reynolds number for y (eqs. (16) or (39)) 

Rej turbulence Reynolds number (eq. (17)) 
j mean rate of strain tensor (eq. (8)) 

S absolute magnitude of mean rate of strain, 7S. 
t time 

mass-averaged velocity vector 
uV mass-averaged fluctuating velocity vector 

u r friction velocity, 

x^ space coordinate vector, i = 1, 2, or 3 
v w surface blowing rate 
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normalized surface blowing rate, — 

i 

closure coefficient in Saffman’s turbulence kinetic energy equation 

closure coefficient in dissipation terms of pseudovorticity or specific 
dissipation rate equation (eq. (14) or (37)) 

closure coefficient in turbulent mixing-energy equation (eq. (13)) 

closure coefficient in decay term of heat-flux equation (eq. (36)) 

closure coefficient in production term of specific dissipation rate 
equation (eq. (14) or (37)) 

closure coefficient in eddy diffusivity expression (eq. (12)) 

Kronecker delta 
eddy diffusivity 

axial rate of strain in homogeneous flow 

e 

turbulence kinetic energy ratio factor at boundary- layer edge, — — 

U e 

Karman constant 

rate of shearing strain in homogeneous flow 

low turbulent Reynolds number factor for closure coefficient in Reynolds 
stress redistribution term (eqs. (16) or (39)) 

closure coefficient in redistribution term of Revnolds-stress equation 
(eq. (35)) 

molecular viscosity 
kiuemat ic viscosity 
mean density 
ins. intaneous density 

inverse of the effective Prandtl number for the turbulent diffusion of 
specific dissipation rate 

inverse of the effective Prandtl number for the turbulent diffusion of 
Reynolds stress or kinetic energy 
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o** inverse of the effective Prandtl number for the turbulent diffusion of 
heat flux 

mass averaged specific Reynolds stress tensor (eqs. (9)) 

T 

r mass averaged specific total stress tensor, molecular plus turbulent 
transport (eq. (6)) 

$ij correlation of pressure and rate of strain fluctuations 
X damping coefficient for \* 

mean vorticity tensor 

uj pseudovorticity or specific dissipation rate 

Subscripts 
w surface value 

o initial value 

e boundary-layer edge value 

Superscripts 

+ "law of the wall" 

+ u + 

u = — , y = y 

UT ’ ; y 


coordinate 

u 

T 

V 
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PROGRESS IN TURBULENCE MODELING FOR COMPLEX FLOW FIELDS 


INCLUDING EFFECTS OF COMPRESSIBILITY 
David C. Wilcox* and Morris W. Rubesin^ 


SUMMARY 


Two second-order-closure turbulence models have been devised that are 
suitable for predicting properties of complex turbulent flow fields in both 
incompressible and compressible fluids. One model is of the "two-equation" 
variety in which closure is accomplished by introducing an eddy viscosity which 
depends on both a turbulent mixing energy and a dissipation rate per unit 
energy, that is, a specific dissipation rate. The other model is a "Reynolds 
stress equation" (RSE) formulation in which all components of the Reynolds 
stress tensor and turbulent heat-flux vector are computed directly and are 
scaled by the specific dissipation rate. Computations based on these models 
are compared with measurements for the following flow fields: (a) low speed, 

high Reynolds number channel flows with plane strain or uniform shear; 

(b) equilibrium turbulent boundary layers with and without pressure gradients 
or effects of compressibility; and (c) flow over a convex surface with and 
without a pressure gradient. The RSE model-computed flow properties generally 
differ by less than 10% from the measurements for all of the applications. 
Discrepancies between two-equation model predictions ana the measured flow 
properties are much larger only for homogeneous flows with suddenly applied 
strain, lending further credence to the notion that a Reynolds stress model 
offers little advantage over a two-equation model for two-dimensional boundary 
layers that are close tc equilibrium. 


INTRODUCTION 


For the past several years, Wilcox and his colleagues have made signifi- 
cant progress in developing second order closure models suitable for computing 
complex turbulent flows, including the effects of compressibility (refs. 1 
and 2) . While early emphasis was focused on two-dimensional attached flows 
with curved streamlines, care was taken to ensure that, at least in principle, 
the models could be applied in a straightforward way to both separated and 
three-dimensional flows. Although most of the modeling effort was devoted to 
the development of a "two-equation" eddy viscosity model of turbulence, this 
work was guided to some extent by the parallel development of a model for the 
Reynolds stress tensor equation. These models were significant improvements 
of the two-equation model originally devised by Saffman (ref. 3) and Wilcox 
(refs. 4 to 6). Most importantly, the model presented by Wilcox and Traci 
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(ref. 6) is as accurate as mixing length theory for equilibrium boundary 
layers* and has a much wider range of applicability, for example, reference 7. 

Throughout the development of these models, the production terms for the 
turbulence kinetic energy were maintained as in Saffman’s original work 
(ref. 3) where the production of turbulence kinetic energy was set proportional 
to the kinetic energy times the absolute value of the rate of strain. Although 
this approach led to a mathematically elegant homogeneous equation for the tur- 
bulence kinetic energy, it introduced modeling for terms that can be derived 
directly. For attached boundary layers, Saffman's additional modeling intro- 
duced no significant differences. When more general flow fields were consid- 
ered, the Saffman-type production terms were demonstrated to contain the poten- 
tial for introducing errors. For example, if the homogeneous turbulence field 
along the centerline of an expanding supersonic nozzle is considered, the 
actual production of kinetic energy is given by 


P 


e 


-(pu-uj) 



( 1 ) 


For this flow field, only normal stresses exist, and thus the stresses that 
appear in equation (1) are positive. As the corresponding rates of strain are 
also positive, equation (1) then indicates a decrease in the kinetic energy 
along the centerline of the nozzle. For this same flow field, the Saffman pro- 
duction term 

P e = a*pSe (2) 

indicates an increase in kinetic energy since the kinetic energy e and the 
absolute magnitude of the rate of strain tensor S are both positive. The 
rather gross discrepancy for this example (and for others as well) suggested 
revision of the production terms in the Wilcox-Traci model. Thus, one objec- 
tive of the current work was to reintroduce the unmodeled production terms 
into the model equation for the kinetic energy. 


During the development of the Wilcox-Traci model, several advances in tur- 
bulence modeling were introduced that are not common, as a group, to most of 
the models found in the literature. First, the effects of compressibility 
have been accounted for through the use of mass-weighted averaged dependent 
variables and some hypotheses regarding terms in which the compressibility 
effects are isolated. Second, several modeling coefficients have been made to 
depend on the turbulence Reynolds number and this, together with the introduc- 
tion of molecular diffusivity as a parallel transport mechanism to the turbu- 
lence, has permitted integration of the transport and modeling equations 
directly to surfaces bordering the flow fields. Thus, the relatively common 
use of a "law of the wall" to define the boundary conditions near surfaces was 
avoided. Finally, a Reynolds stress tensor model for compressible and low 


l A boundary layer is defined as being in equilibrium if its local skin 
friction, momentum and displacement thicknesses are consistent with the 
Ludweig-Tillman formula (ref. 8). 
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Reynolds number flows was developed in parallel with the two-equation model to 
guide the development and assess the limitations of the latter. A secon 
objective of this paper, then, is to modify the production terms xn the Reyn- 
olds stress tensor model and to compare the results of computations based on 
this model with computations from the two-equation model for a series o omo 
geneous and boundary-layer flows that are well documented with measurements of 

the turbulence quantities. 


FORMULATION OF MODEL EQUATIONS 


Mean Flow Equations and Boundary Conditions 

For a compressible fluid in turbulent flow, the conservation equations for 
mass, momentum, and energy are written in terms of the mass-weighted average 
quantities as follows: 


it + ^ ( -“i) ’ 0 

3 .... x . \ = + 3 '~T 


It (pU i } + 3^7 (pu j u i } 


3x ± ■ 317 < P ^> 


it (ph) + ~3 x 7 (pu j h) = n * u i -it + pT Ij 317 - axj (pq j T) 


(3) 

(A) 

(5) 


Here, the symbols tL and qj T denote the specific mass-weighted average total 
shear stress and heat J flux that include the contributions cf both the molecular 
and turbulent transport. These quantities are defined as 
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+ t ij 


( 6 ) 


and 


T v 3h 
q j Pr L 3xj q J 


(7) 


These "specific" quantities are the actual total shear stress and heat flux 
divided by the local mean density. The mean rate of strain tensor appearing 
in equation (6) is given by 


5 ij 


_ 1 (^i 3u j \ 

2 \3xj ZxJ 


( 8 ) 


Finally, tij and qj are the mass-weighted averaged Reynolds stress tensor and 
heat flux vector defined by 
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(9) 


where p is the instantaneous density, < > denotes the time average of the 
enclosed quantity and the superscript " represents the fluctuating part in a 
mass-weighted average formulation. The boundary conditions for equations (3) 
to (5) at surfaces (x 2 = 0) are 

uj = 0 

u 2 = 0 or v w (xj) 

, . / x 3h _ / 9h \ 

i, - M*i> 3 t 2 ‘ 

All flow variables approach free-stream flow conditions in general flow field 
computations. For the special case of two-dimensional boundary layers, the 
boundary condition? at the boundary-layer edge are at x 2 = (Hxj) 

ui = U e ( Xl ) \ 

} (ID 

h = h e (xi) I 

In this section, two models are developed for evaluating the and 

q. defined by equations (9). The first is of the "two-equation" type employ- 
ing an eddy diffusivity that depends algebraically on the dependent variables 
of the two modeling equations. This model is a direct extension of the model 
of reference 6 through the use of unmodeled production terms and a constitu- 
tive relationship that permits nonaligned stress and strain. As with all eddy 
diffusivity models, sudden changes in the rates of strain result in immediate 
changes in the Reynolds stresses and heat fluxes, although the eddy diffusiv- 
ity itself has not had time to respond. To remedy this situation, Saffman 
(ref. 9) developed a method for a two-equation model that permits stress and 
strain to develop at different rates. Although this method showed promise, it 
complicated the modeling process by introducing an additional lag equation. 

The authors felt that the use of Reynolds stress modeling equations could 
achieve these same objectives regarding the relaxation of equilibrium between 
stress and strain in a more direct fashion. The Reynolds stress approach 
eliminated the need for an assumed form of constitutive relationship between 
stress and strain. Further it was felt that an understanding of the rates of 
interaction between the individual Reynolds stresses, not forced by a consti- 
tutive relationship, would lead tc a more general model for three-dimensional 
flows and, perhaps, ould provide guidance for improving the form of the sim- 
pler model and defin the flow conditions for its validity. 
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Two-Equation Model 


The two-equation eddy diffusivity model proposed herein is a direct 
extension of the Wilcox-Traci model of reference 6. The differences occur 
within the production terms of the kinetic energy and specific dissipation 
rate equations 2 and in the constitutive relationship between the Reynolds 
stresses and the mean rates of strain and vorticity. For the reasons cited in 
the Introduction, the production term in the kinetic energy equation is 
retained as the sum of the products of the local Reynolds stress times the 
appropriate mean velocity gradient, equation (1). The analogous term of the 
specific dissipation rate equation expressed in w 2 is modeled as the product 
of the production term of the kinetic energy equation, the ratio w /e and a 
new modeling coefficient. This approach is in keeping with the customary way 
of arriving at the individual terms for the second equation in two-equation 
models. Other than the production terms and the choice of density scaling for 
w, the model equations are virtually identical to those of Wilcox-Traci 
(ref. 6). The definition of the eddy diffusivity is 

c = Z (12) 


where the turbulence kinetic energy and specific energy dissipation rate are 
given by the turbulence modeling equations: 
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where the length scale is represented by 


. 1/2 


l = 


( 13 ) 


The modeling closure coefficients employed are as follows: 


2 These equations differ from what has become the classical two-equation ^ 
model of turbulence (ref. 10) in that the authors interpret e as a "mixing 
energy more akin to 9/4 v' 2 than to the total turbulence kinetic energy in a 
boundary layer, and further, feel that u>, the rate of dissipation of k i net ^ 
energy per unit of kinetic energy, may be a more significant quantity than the 
rate of turbulence dissipation itself. 
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The Reynolds number of turbulence is given by 

B el/2£ 
Re T = 


( 16 ) 


(17) 


The boundary conditions appropriate to these modeling equations when applied 
to boundary layers have been guided by asymptotic analysis and reference to 
other models. The surface boundary conditions for equations (13) and (14) are 
as fellows: 


At x 2 = 0 


or 


where 


e = 0 

20v 

u -*■ — for smooth impervious wall 

6x 2 2 


S U-j 2 

(o •> for a rough and/or porous wall 


/g* v w 
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V (1 + V ) 
w w 


(See ref. 6 for details of the porous or rough wall formulation.) 


At x 2 = 6 


(18) 


e = ’.U e 2 ( Xl ) 
i = 0.09 g* l/4 6(x 1 ) 


(19) 


As discussed in reference 6, the quantity £/6* 1 ^ 4 behaves much like the clas- 
sical mixing length. With this understand ing , the proportionality coefficient 
of 0.09 in equations (19) is readily seen to be consistent with the Escudier 
(ref. 11) eddy-viscosity model for boundary layers. 
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To model Reynolds stresses that do not necessarily align with the mean 
rates of strain, the constitutive relationship relating these quantities was 
expanded to 


ij 


= “3 e6 ij 


+ 2 e 


MS-.) 


9 (8*u) 2 


+ 2S S ) 
mn nnr 


(s. n . + s. a .) 

v im mj jm mx' 


( 20 ) 


where the third term on the right was absent in the Wilcox-Traci model. The 
vorticity tensor used here is defined as 



/ 3u 

\3x. 
X J 


i 



( 21 ) 


The form of equation (20) was guided to some extent by Saffman’s constitu- 
tive relationship, derived in reference 9, where the Reynolds stress is 
expressed in a tensor series to quadratic terms in the mean rates of strain 
and the vorticity 
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( 22 ) 


0) 


where A, B, C, D, F, and G are modeling coefficients. Note that equa- 
tion (22) applies to incompressible flows. In the current work, Saffman's 
distinction between equilibrium and nonequilibrium kinetic energy is dropped, 
and this has consequences in the evaluation of the modeling coef f icients A 
through G. The quantity e in equation (22) is identified here with tne 
kinetic energy of the modeling equation (eq. (13)). 


Equation (22) can be simplified by requiring it to conform with certain 
fundamental experimental observations. The experiment of Tucker and Reynolds 
(ref. 12 ) where a homogeneous, nearly isotropic turbulence field is distorted 
by mean normal strains, indicates that under the normal straining 

uj = exj , uo - -ex ? , and U 3 = 0 (23) 

the normal Reynolds stresses are related approximately as 

f 3 3 = y (in + t ? ?) (24) 


If the conditions represented by equations (23) and (24) are introduced into 
equation (22), it is found that C = 0 is necessary. In addition, in refer- 
ence 13, it was found that a field of homogeneous turbulence in rigid body 
rotation decays without developing an isotropy. This observation requires 


G = 0. 
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If equation (22), with C = 0 and G = 0, is applied to a shear layer 
where the only velocity gradient that exists is 3 ui/3x 2, it is found that 
equation (24) again applies and that B + F = 0 is required to assure that 
the magnitude of the trace of the normal Reynolds stresses is equal to twice 
the kinetic energy. With these values for the coefficients, the Saffman for- 
mulation equation (22) reduces to 


3-3 


- efi . . 
3 ij 
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It should be noted that in a shear layer or in a boundary layer, the last tern 
of equation (25), containing the factor B, is negligibly small. As the 
authors' interest is primarily in such shear flows, the term containing B 
was dropped in the development of equation (20). Equations (25) and (20) 
differ in other respects. The term that appears in equation (20) was 

necessary in the computations to avoid sign difficulties at the outer edge of 
the boundary layers where the values of the normal Reynolds stresses are rela- 
tively small. The term Su^/Sx^ was introduced to account for compressibil- 
ity and to force the trace of the second term on the right in equation (20) to 
vanish. The coefficient A was identified with the y* in equation (12). 
Finally, the quantities 6* and 8/9 were introduced in the last term to 
satisfy some boundary-layer needs. These can be seen by writing equation (20) 
for the logarithmic region of a flat plate boundary layer in incompressible 
flow. There the only mean velocity gradient of consequence is 3u!/3x 2 so 
that 
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with the rest of the components of the rate of strain tensor being orders of 
magnitude smaller. Similarly, the only components of vorticity of importance 
are 


^12 


1 
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2 9x2 


(27) 


With equations (26) and (2/), t i 2 from equation (20) can be written a? 


112 


3u, 
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^ 9 x 2 
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cu 3x2 


(28) 


In this region of the boundary layer, a balance exists between the production 
and dissipation of turbulence kinetic energy in equation (13), thus 


3u i 

pTl2 3x7 ~ 6 * pwe 


(29) 


which leads directly tc 



( 30 ) 
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vhen y* is set equal to unity because of the high turbulence Reynolds number 
in the logarithmic region of the boundary layer (refer to eqs. (16)). The 
normal components of the Reynolds stress tensor in the logarithmic region 
become 


Til 

T22 
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Also, in this region of a f 1, 
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so that equations (29) and (30) combine to indicate 

T /p 

W 

g*i / 2 


(32) 


(33) 


In parallel to these equations, the Reynolds heat flux vectors are repre- 
sented as 


q i 


e 3h _ y* £ 3h 
Pr^, 9x Pr^. 10 3xj 


(34) 


where Fr^, the turbulent Prandtl number, is usually assigned a value of about 
0.9. Comparisons of these modeling conditions with several sets of boundary- 
layer data are given in later sections of this report. 
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Reynolds-Stress Model 

The Reynolds stress formulation presented herein is a straightforward 
extension of the model developed in reference 2. For ' mpleteness, discussion 
in this section reviews some of the material presented in reference 2, partic- 
ularly that pertaining to closure approximations. The key theoretical advan- 
tages for using the Reynolds stress model over the usual eddy viscosity model 
lie in removing the postulates that the principal axis of the Reynolds stress 
tensor align with that of the mean strain-rate tensor and chat sudden changes 
in the mean strain tensor are reflected immediately in the Reynolds stress. 

In terms of mass-averaged dependent variables, the Reynolds stress equation or 
RSE model is expressed in terms of components of Reynolds stress as 
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The components of the Reynolds heat flux are modeled with 
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The specific dissipation rate is given, again, by 
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The modeling coefficients are: 
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Again, the Reynolds number of turbulence, Re x , is given by equations (15) 

and (± 7 ) . At a solid surface, equations (18) again apply and, in addition, 
at x 2 = 0 


q j =° 


T ij 


= 0 


(40) 


At a boundary-layer edge, in addition to equations (19), it is required that 
at x 2 - 6 


^ = 2 I 

T ij = y l ^e 2 ( x l)^ij J (41) 

where the latter assumes isotropic free-stream tui ulence. 

To explain the closure approximations employed, attention is first focused 
on the heat flux equation (eq. (36)), and the specific dissipation-rate equa- 
tion (eq. (37)). The first two terms on the right side of equation (36) are 
the exact heat flux production terms. The next term is the modeled dissipa- 
tion, while the last term is the modeled turbulent diffusion. It is assumed 
that the dissipa:ion of the heat flux is proportional to the heat flux vector 
and that the turbulent diffusion of the heat flux is a simple gradient diffu- 
sion process. For the specific dissipation rate equation, the only difference 
between equation (37) and that used for the two-equation model (eq. (14)), is 
in the form of the coefficient y. An explanation of these differences is 
given in the section entitled "Near Wall Viscous Modifications." 


In the Reynolds stress equations, the first two terms on the right-hand 
side are the exact, unmodeled production terms. The third term is the modeled 
dissipation. The next three terms model the correlation of the pressure and 
rate of strain fluctuations. The final term is the sum of the exact molecular 
diffusion and modeled turbulent diffusion of the Reynolds stress. 


It is illuminating to compare those modeled terms with other models. For 
incompressible flow, the correlation of fluctuations in the pressure and rate 
strain, namely. 
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has been modeled by Launder, Reese, and Rcdi (ref. 14) with the rather complete 
form 


'ij 


C j 0*io yt 


ij 


V 


) ■ <K f 


ij 


P6 


ij 


) - e ( [ 


ij 


! ps i 


j)- 


yeS 


ij 


(43) 


where C t , a, 8, and y 


are the Launder et al. closure coefficients, and 
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(44) 
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The first of the terms on the right of equation (43) is the Rotta “ten- 
dency toward isotropy” term (ref. 15), involving the departure from isotropy of 
the normal Reynolds stresses. The remaining terms represent the interaction 
between the Reynolds stresses and the mean flow. While some of the elements 
in these terms (refer to equations (44)), are similar to the unmodeled produc- 
tion terms, the terms are structured to vanish when the trace of is 

taken, as would be required by continuity for an incompressible fluid within 
equation (42). The relationship of these terms to production and dissipation 
of the Reynolds stresses can be demonstrated most directly by introducing the 
mean rate of strain and vorticity tensors (eqs. (8) and (21)), respectively. 

For example, the exact production terms in equation (35) can then be rewritten 
as 


3u^ 3u^ 
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(45) 


If the strain and vorticity tensors are introduced into equation (43), and 
this general pressure strain correlation is introduced into equation (35), 
there results for incompressible flow 
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This form of the Reynolds stress equations is equivalent to that given in 
reference 16, but the coefficients used in each case have slightly different 
definitions which causes some apparent sign differences. 


From equation (46) it is seen that the pressure rate of strain correlation 
terms, identified oy the modeling constants C lf ut, and A, behave in parallel 
with the production and dissipation terms; the a and 3 complement the factor 
-1 in the first two production terms resulting from strain and rotation, and 
part of can be considered to account for some anisotropic dissipation. 

When equation (46) is applied to a flow with uniform shearing strain, where 


12 


( 47 ) 


1 3u l \ 

Sl2 _ 2 3 ^ ‘ S ?1 

_ 1 3ui 

fll2 ‘ 2 3^7 ' ~ fi 21 

are the only nonzero elements of the strain and vorticity tensors and when 
convection and diffusion are neglected, there results 


and 


Ti * . 2 _ -4 + 4a - 20 

e 3 3Cj 

T22 2 2 - 2a + 40 

e 3 3C X 


X 33 2 2 - 2a - 2g 

e 3 3C! 


The corresponding shearing stress is 


T12 

e 




(n - l)(i 22 /e) + 0(i n /e) + (9/2) 

C[ 


(48) 

(49) 


(50) 


(51) 


Because equations (48) through (50) are not independent, these equations, 
together with equation (51), cannot in themselves be used with specified 
values of the normal and shear Reynolds stresses to define the four constants 
a, 0, y, and Cj . Launder et al. (ref. 14) avoided this difficulty by follow- 
ing Rotta's (ref. 15) suggestion regarding tensor symmetries of the pressure 
strain terms and found that a, 8, Y could each be expressed in terms of a 
pair of constants, Cj and C 2 . They then employed experimental data of the 
normal Reynolds stresses in nearly homogeneous shear flows to define Cj and 
C? through equations equivalent to equations (48) through (50). The approach 
in the present analysis is to require consistency between equations (48) 
through (51) with equations (31), (32), and (33). In ^jdition, it is necessary 
to introduce the assumption that rotation alone does not cause a redistribution 
of the Reynolds stress. This assumption is consistent with the argument that 
eliminated the last term of equation (22). The Reynolds stress model of the 
present analysis, therefore, has been forced to be consistent with the previ- 
ously developed two-equation model, where the normal stresses in the logarith- 
mic region of the boundary layer were set in the ratio Th:t 2 2 : T33 = 4:2:3. 

The values of a, 0, y» and Cj resulting from this model are given in table 1 
along with the values used by Launder et al. (not including special near-wall 
effects) . 
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TABLE 1.- COMPARISON OF CLOSURE COEFFICIENTS IN REYNOLDS STRESS MODELS 


Model 

Coeff J 

Lcient 

C 1 

a 

8 

Y 

Present study 
Launder, Reese, and Rodi 

4.5 

1.5 

0.5 

.76 

0.5 

.11 

1.33 

.36 


The values of the Reynolds stress components in the logarithmic region of 
a flat plate boundary layer utilized in the present model are compared to the 
corresponding values of the model of reference 14 in table 2. 


TABLE 2.- REYNOLDS STRESS COMPONENTS IN THE LOGARITHMIC REGION 
OF A FLAT PLATE BOUNDARY LAYER 


Model 

Reynolds stress component 

T U /e 

T 22 /e 

T 33 /e 

ti 2 /e 

Present study 

Launder, Reese, and Rodi 

-0.89 
- .93 

-0.44 
- .46 

... . 

-0.66 
- .61 

0.31 

.35 


A comparison of these tables reveals that similar values of Reynolds 
stress components can be achieved in a nearly homogeneous shear layer or in the 
logarithmic region of a boundary layer with rather disparate sets of modeling 
coefficients based on different physical arguments. For example, the physical 
assumption that mean flow rotation has no effect on the redistribution of the 
Reynolds stresses force a * 0, and the consequence of this is felt in all the 
other coefficients. One result is that the present model places greater empha- 
sis on the Rotta term, through the larger C lt than does the Launder model. 
Obviously, these models must be tested against data in flows that are differ- 
ent from a nearly homogeneous shear layer to determine which model has the 
better universal character. It is interesting to note that, in terms of the 
symbols used in the present model, Cj = With these modeling constants 

it can be shown that in a uniform shear flow 


2 ^,2 

98* u) 3x 


(52) 


Equation (52) is equivalent to an eddy diffusivity formulation and provides a 
justification for the use of such a concept in nearly uniform shear layers. 

It is interesting, however, that the "eddy diffusivity" is proportional to 
1 2 o rather than to e as in the two-equation model. This difference proves 
extremely important for flows over curved surfaces (ref. 17) and is exploited 
in a later section of this report. 
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APPLICATION 


TO TURBULENT FLOWS WITH PLANE STRAIN AND UNIFORM SHEAR 


In the first series of tests of the two revised models, consideration is 
given to the distortion of high Reynolds number (Re^ >> I) grid turbulence 
caused by an application of plane strain and uniform shear. These cases act 
as tests of the models where near-wall viscous effects are unimportant. For 
these flows the velocity gradient tensor simplifies to 



e k 0 

0 ef 0 

0 0 -(1 + f)E 


(53) 


where c and ic are the axial strain rate and shearing rate, respectively. 
The 11 cases considered are indicated in table 3. 


TABLE 3.- PLANE STRAIN AND UNIFORM SHEAR TEST CASES 


Reference 

l(sec ~ l ) 

<(sec *) 

f 

&*>o 

t max 

Tucker-Reynolds (ref. 12) 

4.45 

0 

.0 

0 

21.5 

0.400 

Reynolds-Tucker (ref. 18) 

4.80 



0 


.220 


14.00 



0 


.135 


12.00 



-1/2 


.145 


3.25 



1 


.255 

Townsend (ref. 19) 

9.41 



0 

70.00 

.150 

Marechal (ref. 20) 

19.00 



0 

60.00 

.135 

Uberoi (ref. 21) 

7.62 



-1/2 

13.5 

.180 


18.20 



-1/2 


.120 


32.50 



-1/2 


.086 

Champagne-Harris-Corrsin 

0.00 

12 

.9 

0 


. 176 

(ref. 22) 








The assumption that the turbulence is homogeneous in a coordinate system 
that moves with the fluid volume, leads to the expression of the two-equation 
model in the following simple form: 


j- - {<-u'v’> ic + [<w ,2 > (1 + f) - (v' 2 > f - <u' 2 )]c} - 8*oue 


(54) 


dw 2 

dt 





e 


{<-u’v’> < 


+ 


[<w' 2 > 


(1 + f) - <v' 2 > f - <u' 2 >]c} - Bu) 3 


(55) 
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where 


and 


<-u t v , > 

e 


jc 

u) 


<u ,2 > 

e 

<v ,2 > 

e 

<w' 2 > 
e 


X 2 


2 l _4_ (k/oj) 2 

3 Z u 98* i + x 2 

2 « £ f 4 (k/cd) 2 

3 u 9B* 1 + x 2 

| + 2 - (1 + f) 

4 ( 7/.)- (1 + f + 1 -) + ( k / co ) 2 


(56) 

(57) 

(58) 

(59) 

(60) 


Note that in writing equations (54) to (60) the turbulent Reynolds number is 
assumed very large so that viscous effects are ignored. The RSE model simpli- 
fies in a similar manner, although the equations are more lengthy; for brevity 
the simplified RSE model equations are not presented. 

The computations of the ordinary differential equations have been per- 
formed with a fourth-order Runge-Kutta integration scheme. To initiate the 
computations, a value for the dissipation rate, co, is needed at t = 0. For 
the Tucker-Reynolds experiment, the initial value of u> has been determined 
from the measured variation of e upstream of the straining region. That is, 
in the absence of strain, e is given by 

= f 1 + \ *Wef ( c ~ c ref)] ( 61 ) 

where subscript ref denotes reference conditions. From the known variation 
of e with t = x/U, the value of ^ re f and hence co at the initial station 
has been determined to be 

Bou 0 = 21.5 sec -1 (62) 

To set 6co 0 for the other flows, we have used values similar to those used by 
Saffman (ref. 9) but with slight differences. 

Of the eleven cases considered, the Tucker-Reynolds (ref. 12), Champagne- 
Harris-Corrsin (ref. 22), and Townsend (ref. 19) experiments are the most com- 
pletely documented. Hence, we first focus on these three cases. 


Tucker-Reynolds Plane Strain Flow 

Figure 1 compares the computed and measured distortion parameter, K, for 
the Tucker-Reynolds plane strain flow; by definition 
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r 

► 


i 


t 


_ <w ,2 > - <u ,2 > (63) 

■ <w' 2 > + <u ,2 > 


As shown, good agreement has been obtained between computed and measured values 
of K when the RSE model is used. Note that the postulated variation of 
with the strain parameter x (see eqs. (39)), viz, 

A* - [f - f exp(-5X)]g* (64) 

is needed to insure the proper relaxation of the anisotropy after the strain 
is removed. Although a less rapid variation of A* with x is acceptab e 
for this flow, equation (64) has been found most suitable overall for all the 
flows considered. Figure 2 shows how close the RSE model predicted normal 
Reynolds stresses are to the corresponding measured stresses when the initially 
anisotropic character of the flow is accounted for. 


Results obtained using the two-equation model are not in as close agree 
ment with the data. At the onset and termination of strain, discontinuous 
jumps Ln the stresses occur. For example, when the strain is removed, the two 
equation model unrealistically predicts an instantaneous return to isotropy. 
Also, the predicted value of K at the end of the straining region is only 
0.39 compared to the corresponding measured value of between 0.55 an 


As a numerical experiment, the two-equation model computation was contin 
ued much farther downstream without removing the strain. The asymptotic va ue 
of the distortion parameter is 0.55, a value nearly identical to that predicted 
by the RSE model. Additionally, it was found that when the Saffman production 
terms are used the asymptotic value is 0.675. With either type of production 
term, the asymptotic value is not achieved until about x 7.6 m. e 
results of this numerical experiment and the fact that the two-equation model 
predicts an abrupt return to isotropy when the strain is removed illustrate 
the fact that even with the new constitutive relation ( eq .. (20)), the two- 
equation eddy viscosity model is limited to near -equilibrium flow conditions. 


Champagne-Harris-Corrsin Uniform Shear Flow 

Figure 3 compares computed and measured Reynolds stresses for the Cham- 
pagne, Harris, and Corrsin uniform shear flow. The entire region of develop- 
ing flow is considered, not merely the asymptotic values. As shown, the RSE 
model predicted stresses that are in relatively close agreement with the data. 
At x = 3.2 m the computed shape parameter, <-u v >//<u - > <v 2 >, of 0.47 is 
•>% lower than the measured value. The only noticeable disagreement occurs 
between computed and measured values for <w’->; the model predicts that as the 
flow approaches equilibrium, 

A 2 (65) 


while the data indicate a ratio of 0.574, somewhat closer to values used by 
Launder, Reese, and Rodi (see table 2). This is a consequence of the fact 
that the model predicts that in equilibrium the normal stresses lie in the 

ratio 
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<u’ 2 >:<v ,2 >:<v’ 2 > = 4:2:3 


( 66 ) 


while the data indicate 


<u ,2 >:<v ,2 >:<w ,2 > = 4:2:2. 4 


(67) 


As with the Tucker -Reynolds computation, the two-equation model predicts a 
slower approach to equilibrium than does the RSE model and the data. A S a “ 1 ’ 
is instructive to examine the asymptotic solution for the two-equation model 
with the conventional and with the Saffman production terms. We find the 
following: 

1 0.48, conventional production 

( 68 ) 

0.76, Saffman production 


and 


<u’ 2 >:<v' 2 > 


:<w / 


1 4:2:3, conventional production 

4: 1.5: 3, Saffman production 


(69) 


Compared with equation (67), equation (69) provides another argument against 
using the Saffman production terms. On the one hand, using the Saffman pro- 
duction terms yields a normal stress ratio of 4:2:3 for a flat plate boundary 
layer (FPEL) . This ratio occurs because the turbulent mixing energy produc- 
tion and dissipation terms balance. However, for the Champagne, Harris, and 
Corrsin flow the balance occurs in the w equation rather than the e equa- 
tion when the Saffman terms are used; consequently the normal stress ratio is 
distorted. On the other hand, using the conventional production terms leads 
to a production-dissipation balance in the e equation both for the Champagne, 
Harris, and Corrsin flow and for flat plate boundary-layer flow. Since the 
RSE model and the Champagne et al. data also indicate that turbulent energy 
production and dissipation balance, using conventional production terms 
appears to yield a more suitable physical representation for uniform shear 

flow. 


Townsend Plane Strain Flow 

Figure 4 compares the computed and measured distortion parameters for 
Townsend's (ref. 19) plane strain flow. Again, reasonably close agreement 
between RSE computed and measured values of K have been obtained. With the 
two-equation model, predicted values of K are about 25% lower than measured 


Other Plane Strain Flows 

Figure 5 presents a summary of computed and measured distortion parameters 
for the ten plane strain flows listed in table 3. Only RSE model predictions 
are shown in the figure. As shown, the predicted values of K at t tmax 
(the time corresponding to the farthest downstream station) lie within 20% of 


18 



corresponding measured values (see table 3), The two-equation model generally 
predicts values of K which are about 70% of measured. The agreement with 
data of long-time asymptotic values is much better, further illustrating the 
fact that the two-equation model fails to accurately describe the rate of 
approach to equilibrium. 


NEAR-WALL VISCOUS MODIFICATIONS 


The conservation and turbulence model equations presented earlier con- 
tained terms involving the molecular viscosity of the fluid in order to permit 
boundary-layer computations with integration through the viscous sublayer and 
no-slip boundary conditions. In this way, use of artificial boundary condi- 
tions such as the "law of the wall" can be avoided ard the uncertainties of 
such boundary conditions far from flat plate boundary-layer conditions can be 
eliminated. For the mean conservation equations it is sufficient to include 
molecular viscosity in the diffusion terms. Similarly, the use of molecular 
diffusion in the various model equations permits imposition of the no-slip 
condition on the Reynolds stress and turbulent heat flux as well as on the 
mean velocity, that is, 


u i = Qi = T ij * 0 at y = 0 


(70) 


However, without further modification it was found by using the sublayer pro- 
gram developed by Wilcox and Traci (ref. 6) that the smooth-wall value of the 
constant in the law of the wall, C, defined by 


C = 


lim 

(u T y/ v )->« 




(71) 


is approximately 7.0 and 7.5 for the two-equation and RSE models, respectively. 
As these values differ from the commonly accepted value of C between 5.0 and 
5.5, further viscous modifications were indicated. Wilcox and Traci (ref. 6) 
have shown that a proper level of C could be achieved by introducing viscous 
damping of the production terms in the e and u> equation. This result can be 
accomplished in the present models by damping the Reynolds shear stress, which 
is implemented in the revised two-equation model by writing 



where 


(72) 


Y* = 1 - (1 - A 2 )exp(-Re T /R e ) (73) 

To allow a different rate of damping of the uj equation production term, we 
also write 

YY* = Y ot [1 - (1 - X 2 )exp(-Re T /R iiJ ) ] (74) 
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where = 10/9 is the high Reynolds number value of y. The closure coeffi- 
cients A, R e and R (jl) are similar to those introduced by Wilcox and Traci 
(ref, 6). 

For the Reynolds stress model, the Reynolds stress itself is damped by 
increasing A* as follows: 


A* = A* [1 - (1 - A 2 )exp(-Re T /R e ) ] -1 


(75) 


while > is modified according to 


X* “ \S [1 “ (1 “ A 2 )exp(-Re T /R w )] 


(76) 


where three closure coefficients A, R e , and R^ again appear. To set the 
value of A, we demand that in a Elasius boundary layer the models predict 
that turbulent fluctuations are damped for Reynolds numbers below the linear- 
stability theory, minimum-critical Reynolds number, Re c = 9><10 4 , Turbulent 
fluctuations are damped provided turbulent energy production t 9u/3y is less 
than dissipation, £*ooe. This condition is satisfied throughout the boundary 
layer up to a plate-length Reynolds number Re^, given by 


Re x - < 


750 


456 


{ A 2 (l - A 2 /3) ’ 


two-equation model 
RSE model 


(77) 


Requiring Re£ 


Re c yields the following values for A: 
' 1 


A ■ 


11 
l 14 


two-equation model 
RSE model 


(78) 


Following Wilcox and Traci, sublayer solutions have been obtained to determine 
loci of values (R e , R (J ) which yield a smooth-wall value for C of 5.5. As 
shown in figure 6, such a locus exists for each model. To select the optimum 
(R e , R w ) pair, the turbulent energy balance between production and dissipation 
has been compared with Laufer's (ref. 23) sublayer data (fig. 7). Closest 
agreement with the data 3 is obtained by using 


R e = 

i ; 

both models 

(79) 

R c " 

I 2 * 

l 3 ; 

two-equation model 
RSE model 

(80) 


3 To further validate the choices for R e and R , a few flat plate 
boundary-layer transition predictions have been made. Transition location and 
width are roughly the same as predicted with the Wilcox-Traci model. Although 
further testing will be needed, it appears that the two revised models probably 
predict transition as well as the Wilcox-Traci model and that transition modi- 

fications devised by Wilcox (ref. 24) apply to both of these models. 
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EQUILIBRIUM BOUNDARY-LAYER APPLICATIONS 


In these applications, the turbulence models are used to compute three 
well documented near-equilibrium boundary layers under incompressible flow 
conditions. The models are then applied to a compressible, flat plate boundary 
layer. Then, an example of a compressible boundary layer is computed to demon- 
strate the effects of relatively strong adverse and favorable pressure gradi- 
ents. These computations are compared with both surface and profile data. 

Incompressible Flows 

To provide definitive tests of the revised models, flows were selected 
(table 4) which include experimental data for all of the Reynolds stress com- 
ponents. Of the three cases, the low Reynolds number of the Andersen adverse- 
pressure gradient case provides the most rigorous test of the models’ ability 
to simulate viscous effects. 


TABLE 4.- BOUNDARf-LAYER EXPERIMENTS USED TO ASSESS THE TURBULENCE MODELS 


Flow 

Data sources 

Flat plate boundary layer 

Karman-Schoenharr (ref. 25) 


skin-friction correlation; 


Klebanoff (ref. 26) data; 


Wieghardt (ref. 27) data. 

Bradshaw adverse pressure gradient 

Bradshaw (ref. 28) data; 


Coles (ref. 27) version of 


Bradshaw data. 

Andersen adverse pressure gradient 

Andersen (ref. 29) data. 


Figures 8 to 10 compare computed and measured shin friction, cp, and shape 
factor, H, along the three equilibrium boundary layers. Over these three 
experiments, both models predict values of skin friction that are generally 
within 5% of the data on both the flat plate and in adverse pressure gradients. 
The predictions are lower than the data, with the RSE model yielding slightly 
lower values than the two-equation model. The computed shape factors are 
within 3% of the measured values for the two-equation model and within 5% when 
the RSE model is used. 

The predicted velocity profiles for both models are compared with the 
data of the three experiments in figures 11 to 13. In figure 11 (the flat 
plate boundary layer), the computations show a slightly thinner boundary layer 
than given by the data, with the larger differences occurring at the low Reyn- 
olds number. Generally, the values of u + are within 5% of the measured 
values For the test cases with adverse pressure gradients, the computations 
yield velocity profiles that are generally less full than the data, but within 
20% of the c nperimental values of u^ even for the poorest comparisons (with 
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the original Bradshaw data or Anderson's 1.8 m station). In these figures 
there is little choice between the two-equation and Reynolds stress models! 

The computed and measured Reynolds stresses are compared in figures 14 
to 16, where the mean velocity profiles are also displayed as u/U 0 versus 

y .?* F °L, al ]t t J iree flows ’ both models predict all four Reynolds stresses to 
within 10% of the measured valves, except in the regions closest to the sur- 
face. This is not surprising when it is recalled that the modeling coeffi- 
cients were based on arguments that in the law of the wall region both the 
turbulent shear and kinetic energy are constant with y (refer to eq. (32) 
and (33) or eq. (51)), whereas the data actually show variations of about 30% 
m e. Changes very close to the surface in the current model, along the lines 
taken by Launder et al. (ref. 14), possibly would improve this situation. It 
is reassuring, however, chat the normal stress component most accurately pre- 
dicted is <v > or T 2 2. and in view of equation (52) explains why the shear 
stress m the present model is predicted rather accurately even though e is 
not. In the normalized mean velocity profiles, the RSE model velocity profiles 
appear to have suffered more deceleration than those of the two-equation model 
On balance, the two-equation model shows slightly closer agreement with the 
boundary-layer data than does the RSE model. 


Compressible Flow 

The turbulence models presented here were extended to compressible flow 
through the assumption that the effects of compressibility can b* introduced 
primarily by use of the Favre mass-weighted dependent variables (ref. 30) in 
both the conservation and modeling equations. Whec Wilcox and Alber (ref. 4) 
first utilized the Favre variables they noted additional terms entered the 
modeling equations that were explicitly dependent on the compressibility of 
the fluid. Also, they suggested that the dissipation rate would require some 
scaling on density to assure proper representation y me compressible terms 
Later, Rubesin (ref. 31) analyzed how the additio ..1 compressibility terms in 
the kinetic energy cou: : be further modeled. Exp.rience with the latter model 
(ret. 32), however, showed that reasonable comparison with data could be 
achieved if these supplemental terms were ignored. This, together with the 
rather successful experience with the Wilcox-Traci two-equation model (ref. 7) 

which also does not use these terms, led to the neglect of extra compressibii-’ 
lty terms in the present model. 

The scaling introduced by compressibility in the present model can be 
illustrated by an analysis of the behavior of the two-equation model in the 
near-wall region of compressible boundary layer on a flat plate. Emphasis will 
be placet, on the iaw-of-the-wall" region, where it is permissible to neglect 
convection and molecular diffusion. The momentum equation (eq. (4)) reduces 
to 


and the model equations simplify to 
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0 ■ «&) 2 


- S‘pue + P* f-y (pe If) 




In the logarithmic region of the boundary layer, the turbulence Reynolds 
number is sufficiently high that 

y* = i \ 

10 } 

Y “ 9 ) 


Also, in this region (ref. 33), 


— = i «.n + B 

u T k 


where 




In this boundary-layer region, a balance exists between the production an is- 
sipation of the turbulence kinetic energy so that the diffusion term in equa- 
tion (82) can also be neglected. From equations (12), (82), (85), and (8b), 
it follows that 


u x 


From equation (81), vith equation (87), it is found that 

{?-)• ' & 
w 

Equations (87) and (88), therefore, represent the density scaling inherent in 
the momentum and turbulence mixing energy model equations. The scaling ques- 
tion would be closed if the forms in the above equations would als satisfy 
the specific dissipation rate equation. When equations (85) through (88) are 
substituted into equation (83) there results 


0 - 


+ 20**1 + 2 o < 


[M^Mv - • •] (99) 


Thus in the limit of the infinitely high Reynolds number where u T> u ? *0, the 
last term vanishes and the usual incompressible t low interrelationship between 
the modeling coefficients results. Therefore, equation (83) is satisfied for 
high Reynolds number, indicating that the density scaling represented by 
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equation in the logiritlmic i re8ion b of h - 'bound^ 11 ? 8 equations and the momentum 
that in this boundfry-layer SSS t^gTa^U ** U 


l = 


e l/2 

e >1/4 


CO 


= 6 * ] 


<y 


( 90 ) 


and identical to that which pyiqi-c a* 4 

wnicn exists in an incompressible flow. 

Plete'a^Maci/numbers'up foVcSle^n” °\? fl *‘ 

number equal to five are eo^ar^ "uh resuUs from tSf “n”? “ a 
in figures 17 and 18 . These formulas have heen *> th ? T Driest 11 formulas 

the bulk of existing data under these conditions to’"bout HO?’ t °J epresent 
between the methods is excellent h..t- h. 4 . ,• out ±10/ * The agreement 

similarity of the density scaling inlir.rlt * Sur * risi *& in view of the 
that which is inherent in the van Driest- f ^ equations ( 87 ) and ( 88 ) and 
largely due to the density scalintr diff* ormul f s * What differences exist are 
the boundary layer. For these flat nlat-I 6 * 1 ^^ . ;Laherent in the w ^ke region of 
and RSE models yield essentially the same resets?" 5 ' the tu0 ~ equatl °“ 

the boundaty-layer ‘ ° T -"^essible flow is 

boundar^layer 1 ^ 'the adiabatic ^r^S “ 

The distribution of the surface nrooc > , y means of a sha P e d centerbody. 
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density nrnf i Ipc q# . _ . ^ provide the mean velocity and 
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Turbulent shear profiles throueh Jhe h^ J he8e 1 data cora P ared to equation ( 85 ). 
shear profiles evaluated Jin; es' T hTs Irlll 5 “?^ 

momentum integrated from the surface to a point within ! a the i mean 
The error introduced by these methods n pn»L * lthi the bou ndary layer, 
models to about ± 15 %. ^ assessment of the turbulence 
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better prediction of the rise in skin friction in the adverse pressure gradient 
region; in the following favorable pressure gradient region, however, they show 
somewhat too large a drop in the skin friction. The two second-order models 
yield essentially equivalent results. 

Figure 21 shows the distribution of the boundary-layer shape factor, 

H = 6*/0, along the test zone. The second-order models yield results that are 
close to the data. The results are generally better than those of the mixing 
length model except for the region immediately downstream of the station where 
the computations were initiated by matching the calculated and experimental 
momentum thicknesses. 

Figures 22 through 27 show comparisons of the calculations based on the 
different turbulence models with mean velocity and turbulent shear profile 
data at the three stations designated in figure 19. Station 1 represents 
essentially flat plate conditions; station 2 is near the end of a long run of 
adverse pressure gradient; and station 3 is at the end of the favorable pres- 
sure gradient region. All of the models well represent the velocity profile 
at static n 1. If one were to be favored, it would be the mixing length model 
which has been fine-tuned to represent equilibrium boundary layers. Figure 23 
shows the turbulent shear profiles computed and measured at this station. 

Near the wall, the mixing length model overpredicts the shear by about 10%, 
whereas the second-order models are about 1C% lower than the measurements.’ In 
the outer part of the boundary layer all the models yield results that are 
about 152 of the peak shear higher than the data. When consideration is given 
to the possible errors in the data, these comparisons are reasonably good and 
demonstrate that the second-order models generally give results in a zero pres- 
sure gradient as good as the fine-tuned mixing length model. 

Figures 24 and 25 show the mean velocity and shear stress profiles at the 
station at the end of the adverse pressure gradient region. The experimental 
mean velocity data have become fuller in the vicinity of the surface, and the 
thickness of the boundary layer has been reduced relative to the data at the 
upstream station. All three models generally reflect this behavior, with the 
two-equation model showing the near-surface behavior best. Surprisingly, the 
RSE model does not show the near-wall fullness in the velocity profile even 
though it provides the best representation of the three models for the skin 
friction and local shear profile. 

After the run of favorable pressure gradient, figure 26 shows the experi- 
mental mean velocity profile fuller and thicker than at the previously indi- 
cated station. All the models represent the experimental data quite well, 
with the second-order models showing a shade better agreement. The experi- 
mental shear stress profile at this station is now about an order of magnitude 
lower than the shear that existed at the previous station. Again the models 
show this drop in shear, however the second-order models overpredict the drop 
in shear in the inner part of the boundary layer. This is consistent with 
their skin-friction values at this station. 

For this example of the behavior of a compressible boundary layer experi- 
encing both adverse and favorable pressure gradients, it is concluded that the 
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second-order closure models, with modeling coefficients that are independent 
of the pressure gradients, represent the experimental data better than a 
mixing-length model with pressure gradient dependent modeling. This conclu- 
sion was also reached in references 38 and 39, where a broad range of compres- 
sible turbulent boundary-layer experiments were compared with the two-e. uation 
model for both attached and separated flows and with the RSE model for attached 
f lows . 


ADVANCED BOUNDARY-LAYER APPLICATIONS 


The models presented here have been compared with the results of boundary- 
layer experiments designed to investigate special aspects of turbulence model- 
ing. The effect of sudden application of transverse shear on the redistribu- 
tion rate between the individual components of the Reynolds stress tensor was 
studied experimentally in references 40 and 41 and compared with our models in 
reference 42. The results of the comparison of the experimental data and the 
computations based on the models showed small improvement over a simple mixing- 
length model but still indicated a too-rapid rate of response to the transverse 
shear. It appears that the success of the RSE model for normally strained 
flows as indicated in figure 1 did not carry over as well to shearing strains. 

Another application of the models has been to the problem of the effect 
of longitudinal curvature on turbulent boundary-layer behavior, the topic 
first recognized by Bradshaw (ref. 43). The Reynolds-stress equations were 
applied to this problem directly through conversion of the coordinates from 
Cartesian to curvilinear with one axis tangent to the body, s, and the other 
normal to the surface, n. The two-equation model, however, required a basic 
reinterpretation of the meaning of the symbol e, the kinetic energy term used 
earlier. 


For flow over a strearawise curved surface with local radius of curva- 
ture R (convex, R > 0; concave, R < 0) the s,n curvilinear coordinate 
system introduces centrifugal and coriolis acceleration terms in the instanta- 
neous equations of motion. In turn, these acceleration terms give rise to 
additional terms inversely proportional to R in the Reynolds stress equa- 
tions. For a thin incompressible boundary layer these equations can be written 

as 
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In these equations, terms with common factors have not been combi led in 
order to reveal the physical source of the individual terms, that is, produc- 
tion, redistribution, etc. The sum of equations (91) through (93) yields the 
kinetic energy equation 
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The specific dissipation rate equation is 
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Equations (91) through (96) were introduced into a boundary-layer code, 
and the results cbtained from this code for comparison with some experimental 
data are given in the following figures. Before discussing these figures, how- 
ever, it is interesting to examine the behavior of the Reynolds stress equa- 
tions in the logarithmic region of the boundary layer where convection and dif- 
fusion terms can be neglected. Note that for this examination centrifugal 
forces are retained. Under these conditions, 
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For a negative value of R, a concave surface, equations (97) to (100) show 
that, relative to the kinetic energy, u* is decreased, v* is increased, w 1 
remains unchanged, and the shear stress t is increased. A convex surface 
results in opposite behavior. 


When equations (95) and (96) were utilized as the model equations for a 
two-equation model in the boundary-layer equations appropriate to longitudi- 
nally curved surfaces, it was found that the additional terms containing R 
did not produce the known curvature effects. This deficiency was corrected by 
exploiting the observations of the similar behavior of the shear stress and 
v 1 in the Reynolds stress equation in order to reinterpret the meaning of e 
for the two-equation model. As in references 1 and 17, e is thought of as a 
’’mixing energy” rather than a ’’kinetic energy” and is defined as 
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With this definition for e and the forms of equations (92) and (95) as guid- 
ance, the e equation for use with the two-equation model only is written, in 
an ad hoc manner, as 
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with 



where e is given by equation (12) and the specific dissipation rate equation 
is equation (96). All the modeling coefficients and relationships employed in 
the two-equation model used earlier are retained. The third term on the left 
side of equation (102) represents the principal extra rate of strain intro- 
duced by the longitudinal curvature. 


These model corrections for the effect of longitudinal surface curvature 
were tested against two boundary-layer flows experimentally measured by So and 
Mellor (ref. 44) over a convex surface. One flow was in a zero pressure gra- 
dient, whereas the second boundary layer experienced an adverse pressure gra- 
dient. Figures 28 through 33 compare the measured data with the computed 
results. Figures 28 and 29 compare measured skin friction and form factors 
with the computed results with and without accounting for the effect of longi- 
tudinal surface curvature. The computed results were matched to the first sta- 
tion by assuming the flow upstream of that station to be on a flat plate of a 
length to yield the correct skin friction there. For the case of the convex 
wall with constant pressure, the unmodified models show none of tne drop in 
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skin friction experienced in the experiment. Except for their behavior imme- 
diately after the first station where the radius of curvature is suddenly 
introduced, the modified models give an excellent representation of the skin- 
friction behavior on the convex surface. This is rather remarkable for the 
two-equation model, when its ad hoc formulation is considered. Both modified 
models also well represent the rise in the form factor experienced by the data. 
For the case of the adverse pressure gradient (fig. 29) the unmodified models 
again do not show an adequate drop in the skin friction, whereas the modified 
models represent the skin— friction data quite well. The models with surface 
curvature corrections also represent the shape factor data in the adverse pres- 
sure gradient well. 

In figures 30 and 31, the computations are compared with the mean stream- 
wise velocity data in law-of-the-wall coordinates. The data correspond to the 
stations farthest downstream where skin friction was measured for the zero and 
adverse pressure gradient cases, respectively. The effect of the convex sur- 
face is to enlarge the contribution of the wake portion of the boundary layer. 
Both the RSE and two-equation models that account for the longitudinal curva- 
ture capture this behavior rather well for both cases of surface pressure dis 
tribution. Again, the added rate of strain in the energy equation of the two- 
equation model yields remarkably close agreement with the data — even somewhat 
better than the RSE model. Note that for the adverse pressure gradient flow 
the large RSE model differences from measurements are mainly due to the models' 
overpredicting c f at the station chosen for comparing theory and experiment. 
Figures 32 and 33 p-ovide a measure of how well each model predicts the various 
Reynolds stresses. Generally speaking, the RSE model predicted stresses are 
as far above the measured values as the two-equation model predicted stresses 
are below measured values. Interestingly, while the near-wall values are more 
closely simulated with the RSE model, the shapes are more closely simulated 
with the two-equation model. 


CONCLUDING REMARKS 


The most significant result of the study is the observation that the two- 
equation model is as accurate ' the RSE model for the class of boundary layers 
considered. Furthermore the new two-equation model constitutive relation per- 
mits accurate prediction of all components of the Reynolds stress tensor. 

The RSE model holds a distinct advantage over the two-equation model for 
turbulent flows which experience sudden changes in the strain rate and/or 
shear. While the two-equation model accurately predicts equilibrium states, 
the rate of approach to equilibrium often is inaccurate. With the RSE model, 
both the rate of approach to equilibrium and the equilibrium state are accu- 
rately predicted. 

Results obtained suggest two tentative conclusions. First, there appears 
to be little advantage in using the RSE model rather than the two-equation 
model for two-dimensional attached boundary layers. Second, for separating 
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flows and for flow subjected to abrupt changes in strain rate and/or shear, 
the RSE model may prove superior to the two-equation model. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, October 15, 1979 
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Figure 1.- Comparison of computed and measured distortion parameter for the 

Tucker-Reynolds plane strain flow. 
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Figure 2.- Comparison of computed and measured normal Reynolds stresses for 

the Tucker-Reynolds plane strain flow. 
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Figure 3.- Comparison of computed and measured Reynolds stresses 
Champagne-Harris-Corrsin uniform shear flow, u - <u'-> / , v ■ 

w = < w ,2 > l /' 
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Figure 4.- Comparison of computed and measured distortion parameter for the 

Townsend plane strain flow. 
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Figure 5.- Summary of RSE model computed distortion parame 

plane strain. 
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Figure 8.- Comparison of computed and measured skin •'riction and shape factor 

for flat-plate boundary- laver flow. 






Figure 9.- Comparison of computed and measured skin friction and shape factor 
for the Bradshaw adverse pressure-gradient flow. 


43 











10 2 10 3 10 4 

u r y/i' 


Figure 12.- Comparison of computed and measured velocity profiles for the 
Bradshaw adverse pressure-gradient flow, x * 2.1 m. 





Figure 13.- Comparison of computed and measured velocity profiles for the 
Andersen adverse pressure-gradient flow. 
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Figure 16. Reynolds stress profiles for the Andersen adverse pressure-gradient 

flow, x = 1.83 m. 
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Figure 17.- Effect of compressibility on the skin friction of an insulated 

flat plate. Re * 10 . 





Figure 18.- Effect of surface cooling on the skin friction on a flat plate at 
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Figure 20*- Skin-friction coefficient distribution in experiment of Levis 

Gran, and Kubota. 
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Figure 22.- Velocity profile at station 1 in experiment of Lewis, Gran 

and Kubota. 
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Figure 23.- Turbulent shear stress profile at station 1 in experiment by 

Lewis, Gran, and Kubota. 
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Figure 25.- Turbulent shear stress profile at station 2 in experiment by 

Lewis, Gran, and Kubota. 
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Figure 26.- Velocity profile at station 3 in experiment by Lewis, Gran 

and Kubota. 
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Figure 27.- Turbulent shear stress profile at station 3 in experiment by 

Lewis, Gran, and Kubota. 
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Figure 28.- Comparison of computed and measured skin friction and shape factor 
for flow over a convex wall with constant pressure. 






Figure 29.- Comparison of computed and measured skin friction and shape factor 
for flow over a convex wall with adverse pressure gradient. 
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Figure 31.- Comparison of computed 
over a convex wall with 
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